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Abstract 


Effective interactions that fit the low energy p-^He experimental data have 
been constructed. They are based on the Resonating Group Method and 
a modified Orthogonality Condition Model in which Pauli and partly Pauli 
forbidden states are removed from the spectrum. Partial waves up to L = 
3 have been considered. The LS force acting between the proton and ^He 
has been included phenomenologically, while the Coulomb interaction has 
been incorporated using a renormalization technique for a screened Coulomb 
interaction. The potentials are also given in a separable momentum space 
form, obtained using the Ernst-Shakin-Thaler (EST) method. In all cases 
the potentials generate phase shifts that fit well the low energy experimental 
data. 

PACS Numbers: 21.45+v, 21.30.-x, 21.60.Gx, 25.10+s 


I. INTRODUCTION 

^He plays an important role in a variety of light-cluster reactions, especially at low 
energies. One such reaction is ^He(d,p)^He which is of interest in primordial nucleosynthesis 
studies and in studies related to the abundance of elements in the Universe. It is also of 
interest as a fusion reaction, as it generates no neutrons. At medium energies, this reaction 
with polarized spins has been used to constrain the deuteron D-state probability |l|. Another 
interesting aspect of this reaction is that the ^He nucleus behaves like a neutron, because 
the proton pair in ^He almost forms a closed shell structure and therefore ^He can serve as 
a good substitute for a neutron target. 

The ^He(d,p)^He reaction can be modeled as a three-body p-n-^He problem which can 
be handled using the Faddeev formalism This presupposes a knowledge of the p-^He 

and the n-^He interactions and of the nuclear potential. The nucleon-trinucleon interaction 
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has been the subject of several investigations in the past. Despite the effort made, the 
interaction is far from being understood. Many questions concerning the range and shape of 
the potential, as well as its energy-, parity-, and ^-dependence. For example, in the optical 
model approach the shape is either pre-chosen (Woods-Saxon p, Gaussian etc.) or is 
htted to specihc scattering cross sections Ba. or to phase shifts which are used as input 
to construct phase equivalent local potentials via inverse scattering methods 0. Some 
experimental data on p-^He elastic scattering up to ~ 10 MeV were analyzed within the 
framework of a separable potential model, phase shifts for the p-^He and n-^H systems were 
generated [pUd . 

The quality of scattering data used as input and their reproducibility is another source of 
ambiguity for the held practitioner. Extensive phase shift analyses for the nucleon-trinucleon 
system were performed in the 1960s and 1970s |Tl|-jT8|. The most recent was carried out 
by Yoshino et al. [j^ in the Ep = 4.0-19.48 MeV energy region. As was pointed out, in 
Ref. 0, performing phase shift analyses and constructing a potential in a specihc channel 
requires extreme care to insure that the interaction has the correct energy-, parity-, and 
^dependence. The mere use of cross sections is not sufficient to uniquely extract all of these 
dependencies. An alternative scheme is to generate the phase shifts from reliable theoretical 
models, such as the Resonating Group Model (RGM). Within the RGM formalism, Reich- 
stein et al. pO extracted the intercluster phase shifts for this system while Furutani et al. 


2T[| calculated the phase shifts by using the Generator Goordinate Method (GGM). Their 


results are in good agreement with the differential cross section and the polarization data in 
the low energy region. Although the calculated energy spectrum of the T = 1 resonances hts 
equally well the experimental data |]^, these experimental data have been replaced in the 
recent article by Tilley et al. p^]. Finally, we mention that an optical potential approach 
has been employed by Teshigawara et al. to study the n-^He system within the Gomplex 
Scaled Siegert Method. 

In this paper, we employ the p-^He nonlocal RGM interaction and its variant 

version, the so-called Orthogonality Gondition Model (OGM) in which the Pauli forbidden 
states (PFS) are removed p5,p^|. Furthermore, to reproduce the experimental phase shifts 
sufficiently well and to treat the degenerate states stemming from the absence of LS forces 
from the RGM kernel, we have introduced phenomenological LS forces using a technique 


similar to that of Ref. p^. Tensor forces could also be considered for the ^Si-^Di, ^Pi-^Pi, 
^P 2 -^F 2 , and ^D 2 -^D 2 channels. However, the corresponding experimental mixing parameters 
were found to be negligibly small |]^ and therefore they have been omitted. 

At higher energies, the absorptive part of the phase shifts could be important, in which 
case an imaginary part should be included in the potential. However, in a recent phase shift 
analysis of Yoshino et al. [O, neither the reflection parameters for the coupling effects in 


different particle channels nor the imaginary part of the phase shifts above the break-up 
threshold were found to be signihcant. Therefore, we will not consider the construction of 
complex p-^He potentials and concentrate instead on real potentials that ht the low energy 
experimental data well. 

Finally, we have endeavored to represent the potentials, that reproduce the phase shifts 
without Goulomb effects, in a separable form with appropriate form factors and parameters 
by using the Ernst-Shakin-Thaler (EST) method [^. Such potentials are useful in solving 


the three-body Faddeev equations in momentum space. 
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In Sec. II we describe our methods to construct the p-^He potential and in Sec. Ill we 
present our results. In Sec. IV we present the separable expansions for the p-^He potentials 
while in Sec. V we present our discussions and conclusions. Finally, some technical details 
and formulas are shifted in Appendices A and B. 


II. THE p-^He POTENTIAL 


The construction of our p-^He effective potential is based on the RGM formalism, which 
is given in conhguration space in the paper by Reichstein et al. 

RGM potential is 


20|. In operator form the 




( 2 . 1 ) 


where W is the energy independent local and nonlocal potential, E is the relative energy 
between the proton and ^He while M is the norm (-integral) kernel (see Appendix A for 
details). 

A common problem encountered in the RGM formalism is the existence of the PFS which 
can be removed via the OGM technique (see Appendix A). The OGM potential is given by 


yOCM 




{ho + W) 




hn 


( 2 . 2 ) 


where ho is the kinetic energy operator between the proton and ^He. It should be noted 
that is an energy independent potential. 

The above operator equations can be brought into a more concrete form by using the 
eigen-values 7 ^ and the eigen-functions Ui of the norm kernel. 


ij=l V-*- 


li 


--{hoij -|- Wi- 


VI - 7i 


hoij] {Uj 


(2.3) 


with 


AA|G,) =7,|G,) 


(2.4) 


and with matrix elements Wij = {Ui\W \ Uj) and hop = {Ui\ho\Uj). We solved Eq. (|2.4|) to 
obtain the eigen-value 7 ^ and the eigen-function Ui for each state of the p-^He system up 
to a large number i = Ng, beyond which the eigen-functions Ui become highly oscillatory. 
These oscillations make the numerical calculations cumbersome and unstable. 

The eigen-functions Ui can be expressed in terms of harmonic oscillator functions Ui (see 
Eq. (|A18|) ). Since the potential becomes inhnite for a PFS( 7 i = 1) the corresponding 

wave function is removed from the eigen-functions of the Schrodinger equation. For the 
partly PFS (0 < 7 ^ < 1) we note that the eigen-values 7 ^ converge quickly to zero as i 
increases and thus for a sufficiently large number Ns, 


R 


OGM 


E IG)»'«{7i 

ij=Ns 


(2.5) 
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Therefore, by using W = \ Ui)Wij{Uj\ = = 0), we obtain instead of Eq. (p.3|) 

the modified OCM (MOCM) potential 


ij=l ij=Ns 

Ns 


ij=l 7* 


(hop + Wi 




Ns 


= = 0 ) + E \Ui) 

ij=l 


VI - 7i 
1 


hop] {Uj 


VI -7* 


(hop + Wi 


ij) 


VI - 7i 


(hop + hEp) (VjI. ( 2 . 6 ) 


If the number of PFS is Np then from the orthogonality condition between (17j| and the 
physical wave function V), given by Eq. ( |A11|) , we obtain 


Np 

Y\Ui) 

p=i 


Vl -7* 


(hop + hEp) 


VI - 7i 


~ (hop + fEp) 


{U, 


= 0 . 


(2.7) 


Thus, the part of the potential in Eq. (|2.6| ) has no effect in the Schrodinger equation. Con¬ 
sequently, we introduce our MOCM potential as follows. 


Ns 


1/MOCM ^ yRGM(^^ = 0) + ^ \Ui) 

ij=Np-\-\ 


VI -T 


(hop + fEp) 


VI - 7i 


(hop + hEp) (hjl- 

( 2 . 8 ) 


For the £ = 0 state of the p-^He system, the eigenvalues y* are 71 = 0.98993, 72 = 0.12020, 
73 = 0.01373, and so on, i.e. there are no PFS in the system. For higher partial waves, the 
absolute values of 7 ^ are much less than one, and therefore Np = 0 as well. Thus, Eq. (p. 8 |) 
is not only free from a diverging term caused by 7 ^ = 1 , but also from numerical errors due 
to presence of oscillator functions with higher frequencies. 

Each term of Eq. (|2.8|) can be easily transformed into its momentum space using Fourier 
transforms (see Appendix B). It is found that the on-shell parts of the OCM scattering 
amplitude obtained with the transformed potentials are almost the same as those of the 
RGM amplitude, although the off-shell parts differ. For the spin-triplet P- and D- waves 
and with = 20 the resulting phase shift differ slightly. However, the differences are small 
enough and do not give rise to any important effects in the parameter htting of these triplet 
channels. 

Within the above formalism, a realistic force with an LS component has rarely been con¬ 
sidered until now [^|. As a consequence, the intercluster p-^He potential of the MOCM, 
Eq. ( p. 8 D , does not include an LS force either. Although in the RGM formalism the NN 


interaction could include an LS and tensor force, this may necessitate some additional an- 
tisymmetrization procedure. However, if we adopt a proper intrinsic cluster function that 
includes the effects of these forces, we could include these effects in the inter-cluster potentials 


by means of a folding procedure. In previous work McIntyre and Haeberli introduced 
a phenomenological LS force to reduce the degeneracy of the L-induced potentials. In the 
present work we use the same technique but with an extended form for the LS forces. The 
LS dependence in the intercluster potentials for the spin-triplet channels is introduced via 
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^OCM(LS) = {C„ + Cii • S + C2(L ■ Sf + C3(L- s)S}i/>''ocm 


(2.9) 


where Co, Ci, (72,and C 3 are parameters adjusted to fit the experimental phase shifts. For 
now on we shall, for simplicity, use from now on the abbreviation OCM instead of MOCM. 
The parameters for the p-^He interactions are shown in Table.|. 
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TABLES 



TABLE I. 

Parameters of the p- 

-3He interactions. 


State 

Co 

Cl 

C2 

Cs 

3P0 3Pi 3P2 


0.0788 

-0.0740 


3Di 3D2 3D3 


0.0446 

0.0616 


3F2 3F3 3F4 


0.0360 

0.0084 
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Tensor forces should be considered, in principle, for the ^Si-^Di, ^Pi-^Pi, ^P 2 -^F 2 , and 
^D 2 -^D 2 coupled channels. However, as mentioned in the introduction, in the phase-shift 

it was found that the mixing parameters for the p-^He 
= -0.2719 ± 0.0097deg, £7 = -0.4001 ± 0.3490 deg, 
4 = -0.0269 ± 0.2510 deg, and = 0.0827 ± 0.0129 deg at Ep = 19.48 MeV for the ^Si- 
^Di, ^Pi-^Pi, ^D 2 -^D 2 , and ^P 2 -^F 2 partial waves, respectively. Therefore, tensor forces were 
omitted. 


analysis by Yoshino et. al. 
interaction are very small; i.e. 


III. PHASE SHIFT CALCULATIONS 


The p-^He phase shifts were calculated using the program “GSE” developed in Refs. 


30| , |3l|| . The GSE method is one generalization of Bateman’s method in which the on-shell 


and half-off-shell properties of the t-matrix are exact when the on-shell momentum is chosen 
to be the Bateman parameter. In order to take into account the Goulomb interaction for the 
p-^He system, Reichstein et al. introduced a kind of Goulomb force in the RGM formalism 
which takes the form of an Error function [^. In momentum representation, the Error 
function can be expanded in the short range region in terms of Gaussian functions while 
in the long range region it can be expressed in terms of a Goulomb potential between the 
^He cluster and the proton. In order to obtain the phase-shift modihcation by the Goulomb 
interaction, we adopt a screened Goulomb potential of the Yukawa type. The difference 
between the pure Goulomb and the screened Goulomb potential in the long range region is 
corrected by means of a renormalization technique 
The screened Goulomb potential is dehned by 


V^ir) = 


ZZ'e 


'^2 


■ exp(—r/R) 


(3,1) 


where R is called the screened Goulomb range parameter which, in our case, has the value 
of R = 810.753 fm for S-wave, R = 1128.677fm for P-wave, R = 1439.062 fm for D-wave, 
and R = 10158.163 fm for F-wave. The partial wave expansion of U^, in momentum space, 
is given by 


VtHk, k') 


— 2 ^' — 


(3.2) 


where Qe{x) is the Legendre function of the second kind. The calculated screened Goulomb 
phase shift using Eq. (|3.2|) is given by 

6 f{k) = ae{k)-(j){k,R) (3.3) 


where (Ti{k) is the Goulomb phase, (j){k, R) is the renormalization phase given by (j){k, R) = 
//(In 2kR — j) + ■■■, rj is the Sommerfeld parameter r] = laZZ'e^/k, p is a reduced mass, and 
7 is the Euler constant 7 = 0.5772156 ■ ■ ■. 

In order to calculate the Goulomb modihed nuclear phase shifts we hrst put the 
total potential with the short range nuclear part and the Goulomb part into partial 

wave form, 

Ve{k, k') = k') + U/(fc, k'). (3.4) 
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Then, the calculated total phase shift is given by 


= 5!^ + 5f{k) = 5l^{k) + at{k) - (j){k, R). (3.5) 

Consequently the Coulomb modihed nuclear phase shift 5f^{k) is 

y«(4) = - 5f(k) 

— + 4>{k^R) — (3.6) 

The genuine Coulomb modihed nuclear phase shifts are obtained via Eq. (p.6|) in the long 
range limit, 

= lim R)] - Mk) 

R-^oo 

= 6i{k) - ae{k), (3.7) 


where 6e{k) is the total phase shift generated by the nuclear plus Coulomb potentials. The 
results, for the ^Sq, ^Si partial waves, are shown in Fig. and compared with the results 


of the phase shift analyses 
channels are well reproduced 
RGM potentials. 


ri|,|T^,|IJ-0 . It is seen that the phase shifts for these 


This conhrms the on-shell equivalence of the OCM and the 
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FIGURES 




FIG. 1. p-^He phase shifts for the ^Sq and ^Si partial waves. The experimental data are: o 
Tombrello [12] □ Me Sherry and Baker [14], A Morales et al. [16], o Miiller et al. [17], + 
Beltramin et al. [10], and • Yoshino et al. [19]. The dashed line denotes the RGM and the solid 
line the RGM+OGM results, respectively. 

In Fig. the phase shifts for the ^Pq, ^Pq, ^Pi, ^P 2 partial waves are shown. The spin 
singlet ^Pq phases for the RGM and the RGM+OGM potentials are obtained without an 
LS force. In contrast, the ^Pq, ^Pi, ^P 2 partial waves are spin triplet, and therefore the LS 
force is taken into account. It is seen that our results are in very good agreement with the 
experimental data, especially for energies below 10 MeV. This demonstrates the importance 
of the LS force in the potential. 
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FIG. 2. p-^He phase shifts for the ^Pi, ^Pq, ^Pi, and ^P 2 partial waves. Symbols for the 
experimental data are the same as in Fig. 1. For the singlet channel the dashed line denotes the 
RGM and the solid line the RGM+OCM results. For the triplet channel the dashed line denotes 
the RGM, the dashed-dotted line the RGM-I-OGM, and the solid lines the RGM-|-OCM-|-LS results, 
respectively. 


In Fig. ^the phase shifts for ^D 2 , ^Di, ^D 2 , and partial waves are shown. The spin- 
singlet ^D 2 partial wave is calcnlated using the RGM-I-OCM potential without inclusion of 
an LS force. Once again, the results are in good agreement with the experimental data of 
Yoshino et al. in the lower energy region. It should be noted that, unlike for the P wave 
results, the LS force is less important for the ^Di, ^D 2 , and ^Da partial waves. 






Ecm. [MeV] E, „ [MeV] 

FIG. 3. p-^He phase shifts for the ^D 2 , ^Di, ^D 2 , and ^Da partial waves. The notation is the 
same as in Fig. 1 and 2. 


In Fig. the ^Fa, '^F 2 , '^Fa, and '^F 4 partial wave phase shifts are shown. Again, the 
^Fa is a spin singlet. The results were obtained using the RGM and RGM-I-OGM potentials 
only. In order to obtain the relevant constants for the LS force in the ^F 2 , ^Fa, and ^F 4 


channels, we employed the latest experimental data by Yoshino et al. ro 
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FIG. 4. p-^He phase shifts for the ^F 2 , ^Fa, ^F 4 partial waves. The notation is the same 
as in Fig. 1 and 2. 


IV. SEPARABLE EXPANSION OF THE p-^He POTENTIALS 


In order to be useful the potential yOCM^LS) j^iust be converted into a separable expansion 
form. To achieve this we employed the EST separable expansion method |^] which is briefly 
described below. 

The EST rank N separable potential {p,p') is dehned in terms of the form factor 
giip) and the coupling constant \ij 


N 

V^^^{p,p') = ^ gi{p) Xij gj{p). 

The form factor is given by the R-matrix, 

g,{p) = {p\V\^{E,)) = R{p,ki;E,), 


(4.1) 


(4.2) 


where Ei is a hxed energy point and ki is the on-shell momentum at this energy. The 
expansion energies used are Ei=lMeV for the rank-1 case, and Ei=lMeV, E2=10MeV, and 
£'3=15MeV for the rank-3 case. 

The R-matrix satishes the equation 


foo 

R{p,kpE,) = V{p,ki)+V / R(p,p')Go(p;E,)R(p',fcpE,)f-dp' 

Jo 


= V{p^k,)+V [ 

Jo 


~ V{p,p')R{p',ki]Ei) p'^ 

Ei - p'V2/i 27r2 


dp' 


(4.3) 
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where V denotes the principal value integral. For practical reasons the form factor gi{p) is 
expanded in terms of polynomials, 

20 

9nip) = e-^P"p^ Y. , (4.4) 

m=l 

where (3 = 1/0.825 and an,m are htting parameters. 

The coupling constant A^- is dehned by 


rco 

{X-% = {^iE,)\V\^{E,)) = R{h,k,-Ei)+V / Riki,p-,E,)Goip;E.^R{p,kp,E,)^dp. 

Jo 


(4.5) 


A. Singlet channel 


The singlet p-^He potential does not contain an LS force. We endeavored to construct 
two different p-^He potentials, one of rank-1 and the other of rank-3. The rank-1 parameters 
for the form factor are the same as those of the hrst term of the rank-3 but the coupling 
constant is different. 

The parameters for these potentials for the ^Sq state are given in Tables II and III while 
the phase shift results are presented in Fig. It is seen that the difference between the 
phase shifts, especially in the lower energy region, is very small. 


TABLE II. Parameters an,m of the form factor of Eq. 
column lists the rank-1 parameters with An = 2.976586. 

(4.4) for the ISq potential. The first 


^2,m 


0.81290761 xlO^ 

-0.37821814x102 

-0.13596682x102 

-0.13818233x102 

0.45121340x102 

0.14012550x102 

0.26993988x102 

-0.86975431x102 

-0.27226075x102 

-0.29584805x102 

0.78262662x102 

0.21711608x102 

0.22158774x102 

-0.47493102x102 

-0.11624436x102 

-0.12179712x102 

0.20523082x102 

0.42622121x101 

0.52046713x10^ 

-0.68468420x101 

-0.11659454x101 

-0.17809791x101 

0.18744208x101 

0.26053262 

0.49289681 

-0.44089876 

-0.54300629x10-1 

-0.11013105 

0.90888217x10-1 

0.11621725x10-1 

0.19710935x10-1 

-0.16177816x10-1 

-0.23952017x10-2 

-0.27963987x10-2 

0.23984064x10-2 

0.41701960x10-3 

0.31072739x10-3 

-0.28510226x10-3 

-0.56522712x10-^ 

-0.26665392x10-^ 

0.26325182x10-^ 

0.57510182x10-3 

0.17360057x10-3 

-0.18369762x10-3 

-0.43035052x10-° 

-0.83650320x10-'^ 

0.94090317x10-1 

0.23179820x10-1 

0.28751338x10-8 

-0.34035371x10-8 

-0.86968112x10-° 

-0.66353945x10-1° 

0.81822792x10-1° 

0.21475347x10-1° 
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0.91778681x10-^2 

-0.57298542x10-^^ 


-0.11672988x10-^^ 

0.74478363x10-^^ 


-0.31250844x10-^2 

0.20237810x10-^^ 


TABLE III. Coupling constants Xij = Xji of the rank~3 ^Sq potential. 

All =-0.12888040x102 
Ai 2 =-0.12924909x102 
Ai 3 = 0.25638703x102 
A 22 =-0.11179525x102 
A 23 = 0.21826395x102 
A 33 =-0.42790529x102 
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Ec.m. [MeV] 

FIG. 5. Phase shifts for the ^Sq partial wave without the Coulomb effects. The crosses denote 
the RGM+OCM results, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


The results for the ^Pi channel are given in Tables IV and V while the phase shifts are 
given in Fig. 6. As can be seen from this figure the rank-1 potential fails to reproduce the 
results beyond ~ 4 MeV. Similar results were obtained for the ^D 2 and channels. These 
results are presented in Tables Vl-lX and plotted in Figs. 7-8. 


TABLE IV. Parameters an,m of the form factor of Eq. (4.4) for the ^Pi potential. The first 
column exhibits the rank-1 parameters with An = —0.29317504. 



02,m 

(^3,m 

-0.94249747 xlO^ 

-0.80751728x10^ 

-0.54943593x101 

0.16850309x102 

0.72033423x10^ 

0.28748548x101 

-0.25259235x10^ 

-0.90782035x10^ 

-0.33647532x101 

0.22361127x102 

0.46110606x10^ 

0.41169528 

-0.14005952x102 

-0.13188315x10^ 

0.60064787 

0.67212849x10^ 

0.16478322 

-0.41395211 

-0.26035140x10^ 

-0.27097734x10-1 

0.86574730x10-1 

0.83231138 

0.41703584x10-1 

0.28862915x10-1 

-0.22037742 

-0.29197374x10-1 

-0.28632910x10-1 

0.47917437x10"^ 

0.11484498x10-1 

0.11370739x10-1 

-0.84391784x10-2 

-0.29953237x10-2 

-0.29062500x10-2 

0.11860729x10-2 

0.55358368x10-3 

0.52610797x10-3 

-0.13104955x10-3 

-0.74563678x10-^ 

-0.69658876x10-1 

0.11203792x10-^ 

0.73865049x10-3 

0.68075554x10-3 

-0.72718918x10-® 

-0.53624860x10-° 

-0.48896135x10-° 

0.34933248x10-'^ 

0.28090702x10-'^ 

0.25397873x10-'^ 

-0.11963972x10-3 

-0.10294665x10-3 

-0.92452666x10-° 

0.27487008x10-^° 

0.24947694x10-1° 

0.22283655x10-1° 

-0.37800739x10-^2 

-0.35784010x10-12 

-0.31823024x10-12 

0.23427534x10-^^ 

0.22924616x10-11 

0.20314053x10-11 


14 










TABLE V. Coupling constants Xij = Xji of the rank-3 ^Pi potential. 

All =-0.86234010x101 
Ai 2 = 0.26099258x102 
Ai3 =-0.24411307x102 
A 22 =-0.88300863x102 
A 23 = 0.85424554x102 
A 33 =-0.83781354x102 
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FIG. 6. Phase shifts for the ^Pi partial wave without Coulomb effects. The crosses denote 
the RGM+OCM result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE VI. Parameters an,m of the form factor of Eq. (4.4) for the ^D 2 potential. The first 
column gives the rank-1 parameters with An = 0.49171217x10^ 



0-2,m 

0-3,m 

0.68105470 

0.19763270x101 

0.16194259x101 

-0.19186565x101 

-0.54138849x101 

-0.45399709x101 

0.27889420x101 

0.67678908x101 

0.50617192x101 

-0.25046661x101 

-0.50530126x101 

-0.33071540x101 

0.15881034x101 

0.26186515x101 

0.14542318x101 

-0.76207628 

-0.10174428x101 

-0.45838231 

0.28729893 

0.31100450 

0.10636485 

-0.86461811x10-1 

-0.76634068x10-1 

-0.17716414x10-1 

0.20858825x10-1 

0.15393708x10-1 

0.18712562x10-2 

-0.40260171x10-2 

-0.25251882x10-2 

-0.43762336x10-^ 

0.61860736x10-3 

0.33681168x10-3 

-0.26015437x10-^ 

-0.75152221x10-^ 

-0.36220702x10-^ 

0.57462271x10-3 

0.71556851x10-3 

0.31045680x10-3 

-0.69663479x10-° 

-0.52783394x10-3 

-0.20899733x10-° 

0.57050021x10-1 

0.29679582x10-1 

0.10841787x10-1 

-0.32922530x10-3 

-0.12425007x10-3 

-0.42224769x10-° 

0.13358706x10-° 

0.37349976x10-1° 

0.11883187x10-1° 

-0.37020130x10-11 

-0.75924107x10-12 

-0.22720832x10-12 

0.65926451x10-13 

0.93178858x10-1^ 

0.26314644x10-1^ 

-0.66816378x10-13 

-0.51997857x10-1° 

-0.13885275x10-1° 

0.28516854x10-11 
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TABLE VII. Coupling constants = \ji for the rank-3 ^D 2 potential. 

All = 0.23985939x10^ 

Ai2 =-0.13537297x10^ 

Ai 3 = 0.66626420x103 
A 22 = 0.74471421x103 
A 23 =-0.34111641x103 
A 33 = 0.13709213x103 
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FIG. 7. Phase shifts for the ^D 2 partial waves without Coulomb effects. The crosses denote 
the RGM+OCM result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 

TABLE VIII. Parameters an,m of the form factor of Eq. (4.4) for the ^Ea potential. The first 
column lists the rank-1 parameters with An =-0.157512965x10^. 


02,m 


03,m 

-0.10309579 

-0.10491693x103 


-0.10971956x103 

0.25614578 

0.21551217x103 


0.19720076x103 

-0.35470074 

-0.26919503x103 


-0.23113579x103 

0.30827988 

0.20398472x103 


0.15842458x103 

-0.18968846 

-0.10780556x103 


-0.73781489 

0.88041294x10"^ 

0.42456486 


0.24470228 

-0.31866032x10-^ 

-0.12903376 


-0.57371456x10-3 

0.91323183x10-2 

0.30698165x10-3 


0.84797207x10-2 

-0.20840227x10-2 

-0.57388079x10-2 


-0.30810984x10-3 

0.37892090x10-3 

0.84146200x10-3 


-0.20824797x10-3 

-0.54767188x10-^ 

-0.96337919x10-^ 


0.63773382x10-^ 

0.62640005x10-3 

0.85660575x10-3 


-0.10605018x10-^ 

-0.56298706x10-3 

-0.58902484x10-° 


0.11988500x10-3 

0.39349866x10-^ 

0.31325432x10-7 


-0.97107918x10-7 

-0.21061797x10-3 

-0.13004858x10-3 


0.56943638x10-3 

0.84350422x10-^° 

0.42992976x10-3° 


-0.23915836x10-° 

-0.24379678x10-^3 

-0.11462659x10-33 


0.69843623x10-33 

0.47882660x10-33 

0.23655180x10-33 


-0.13400069x10-32 

-0.57032237x10-33 

-0.32743658x10-33 


0.15094646x10-34 

0.31009923x10-37 

0.21752827x10-37 


-0.75136464x10-37 
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TABLE IX. Coupling constants Xij = Xji for the rank-3 potential. 

All =-0.11253575x10^ 

Ai 2 = 0.26062364x10^ 

Ai 3 =-0.14514498x10^ 

A 22 =-0.63200447x10^ 

A 23 = 0.35959447x10^ 

A 33 =-0.20743460x10^ 
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E,^[MeV] 

FIG. 8. Phase shifts for partial wave without Coulomb effects. The crosses denote the 
RGM+OCM result, the dashed line the rank-1 result, and the solid line the rank-3 result, respec¬ 
tively. 


B. Triplet Channel 


The parameters for the spin-triplet partial waves ^Si (^Po,^Pi, ^P 2 ), (^Di, ^ 02 ,^ 03 ), 
and (^F 2 , ^F 3 ,^F 4 ), are given in Tables X-XXIX while the phase shifts are given in Figs. 
9-18. In the case of the ^Si partial wave the rank-3 separable potential reproduces quite well 
the RGM-I-OCM phase shift. In the other cases the rank-3 separable potential reproduces 
equally well the RGM-I-OCM-I-LS phase shifts. 

TABLE X. Parameters an,m of the form factor of Eq. (4.4) for the ^Si potential. The first 
column gives the rank-1 parameters with An =-0.25598940x10^ 


®2,m 

®3,m 

0.74418704x10^ 

-0.13197287x103 

-0.16817810x102 

-0.99634394x101 

0.10665018x103 

0.11212182x102 

0.23859818x102 

-0.30740737x103 

-0.35992162x102 

-0.27331253x102 

0.28417959x103 

0.28919328x102 

0.20995060x102 

-0.17596997x103 

-0.15737689x102 

-0.11699246x102 

0.75921030x102 

0.56602713x101 

0.50516284x101 

-0.25132501x102 

-0.14815321x101 

-0.17466475x101 

0.68293763x101 

0.30459200 

0.48893164 

-0.16061242x101 

-0.56883879x10-1 

-0.11058936 

0.33443137 

0.11406159x10-1 

0.20040057x10-1 

-0.60527031x10-1 

-0.23718437x10-2 

-0.28772181x10-2 

0.91323541x10-2 

0.42859938x10-3 

0.32324420x10-3 

-0.11020710x10-2 

-0.60155849x10-^ 

-0.28013348x10-^ 

0.10296990x10-3 

0.62866235x10-3 

0.18394221x10-3 

-0.72504967x10-3 

-0.47999987x10-° 

-0.89281320x10-'^ 

0.37402881x10-° 

0.26272042x10-1' 

0.30874015x10-3 

-0.13612727x10-'^ 

-0.99973788x10-° 

-0.71610131x10-1° 

0.32922428x10-° 

0.25040321x10-1° 

0.99454633x10-12 

-0.47290138x10-11 

-0.37035448x10-12 
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-0.62302608x10 


0.30443574x10 


0.24485301x10 


-14 


-13 


-14 


TABLE XL Coupling constants Xij = Xji for the rank-3 ^Si potential. 

All = 0.25332051x101 
Ai 2 =-0.44177961 
Ai 3 = 0.33554529x101 
A 22 =-0.19258003 
A 23 = 0.11500883x101 
A 33 =-0.67841091x101 
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FIG. 9. Phase shifts for the ^Si partial wave without 
the RGM+OCM result, the dashed line the rank-1 result, 
respectively. 

TABLE XII. Parameters an,m of the form factor of Eq. 
column lists the rank-1 parameters with An =-0.42316727. 

25 30 

Goulomb effects. The crosses denote 
and the solid line the rank-3 result, 

(4.4) for the 3po potential. The first 


02,m 


03,m 

-0.66322212x101 

-0.79077911x101 


-0.59125005x101 

0.10444487x10^ 

0.55415998x101 


0.24492255x101 

-0.17647766x102 

-0.86712354x101 


-0.38446872x101 

0.16035080x102 

0.44944025x101 


0.77187900 

-0.10440381x102 

-0.17971337x101 


0.35090483x10-1 

0.52745557x101 

0.78617775 


0.12686539 

-0.22025734x101 

-0.51578766 


-0.30787262 

0.77506575 

0.31027724 


0.24158041 

-0.22810778 

-0.13698661 


-0.11373402 

0.54968268x10-1 

0.43590322x10-1 


0.36741567x10-1 

-0.10611771x10-1 

-0.10160062x10-1 


-0.85776771x10-2 

0.16122722x10-2 

0.17572875x10-2 


0.14804394x10-2 

-0.18998460x10-3 

-0.22687698x10-3 


-0.19056704x10-3 

0.17117673x10-^ 

0.21835886x10-^ 


0.18287608x10-^ 

-0.11593237x10-3 

-0.15527041x10-3 


-0.12969524x10-3 

0.57638824x10-'^ 

0.80073563x10-'^ 


0.66726978x10-'^ 

-0.20292626x10-3 

-0.28985088x10-3 


-0.24102829x10-3 

0.47658438x10-1° 

0.69528124x10-1° 


0.57703914x10-1° 

-0.66681851x10-12 

-0.98848850x10-12 


-0.81884000x10-12 

0.41875253x10-1^ 

0.62812081x10-1^ 


0.51931988x10-1^ 
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TABLE XIII. Coupling constants Xij = Xji for the rank-3 ^Pq potential. 

All =-0.32152588x102 
Ai 2 = 0.93930066x102 
Ai 3 =-0.90374767x102 
A 22 =-0.29056296x103 
A 23 = 0.28341791x103 
A 33 =-0.27758530x103 
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FIG. 10. Phase shifts for the ^Pq state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XIV. Parameters an,m of the form factor of Eq. (4.4) for the ^Pi potential. The first 
column exhibits the rank-1 parameters with An =-0.15329984 



0-2,m 

0-3,m 

-0.11983147x102 

-0.13550093x102 

-0.90367514x101 

0.17014168x102 

0.73213312x101 

0.24098258x101 

-0.28697894x102 

-0.11569940x102 

-0.40100766x101 

0.25268582x102 

0.43479794x101 

-0.58369026 

-0.16200758x102 

-0.11052072x101 

0.96722200 

0.81469516x101 

0.52139426 

-0.13854134 

-0.34392479x101 

-0.63361772 

-0.38566501 

0.12378490x101 

0.47703657 

0.35442445 

-0.37446386 

-0.22713933 

-0.17258809 

0.92646826x10-1 

0.74443476x10-1 

0.56334788x10-1 

-0.18290727x10-1 

-0.17583617x10-1 

-0.13201035x10-1 

0.28292359x10-2 

0.30615614x10-2 

0.22815526x10-2 

-0.33808968x10-3 

-0.39667201x10-3 

-0.29380436x10-3 

0.30793834x10-^ 

0.38250521x10-^ 

0.28190681x10-^ 

-0.21030263x10-3 

-0.27223565x10-3 

-0.19982189x10-3 

0.10522448x10-3 

0.14041899x10-° 

0.10271462x10-° 

-0.37222405x10-3 

-0.50808201x10-3 

-0.37053918x10-3 

0.87716868x10-1° 

0.12175644x10-° 

0.88550576x10-1° 

-0.12300073x10-11 

-0.17282223x10-11 

-0.12535112x10-11 

0.77324366x10-1^ 

0.10955268x10-13 

0.79236297x10-1^ 
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TABLE XV. Coupling constants Xij = Xji for the rank-3 ^Pi potential. 

All =-0.10202649x102 
Ai 2 = 0.31235498x102 
Ai 3 =-0.33737459x102 
A 22 =-0.10318492x103 
A 23 = 0.11347250x103 
A 33 =-0.12541685x103 
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FIG. 11. Phase shifts for the ^Pi partial wave without Coulomb effects. The crosses denote 
the RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XVI. Parameters an,m of the form factor of Eq. (4.4) for the ^P 2 potential. The first 
column presents the rank-1 parameters with An =-0.66395138x10“^ 



®2,m 

03,m 

-0.18802681x102 

-0.20111944x102 

-0.11866155x102 

0.24695039x102 

0.86908371x10! 

0.19942254x10! 

-0.41611877x102 

-0.13903725x102 

-0.36405995x10! 

0.35688295x102 

0.31275599x10! 

-0.22992102x10! 

-0.22579637x102 

0.31524135 

0.20651004x10! 

0.11307515x102 

-0.49795244x10“! 

-0.47700417 

-0.48182932x10^ 

-0.68502022 

-0.42575885 

0.17672518x101 

0.64872070 

0.44814211 

-0.54652963 

-0.32730654 

-0.22367571 

0.13797319 

0.10952640 

0.73558687x10“! 

-0.27690659x10“! 

-0.26106999x10“! 

-0.17285137x10“! 

0.43380605 xl0“2 

0.45661980x10“2 

0.29907200x10“2 

-0.52342569 xl0“3 

-0.59304311x10“® 

-0.38526921x10“® 

0.48023975x10“^ 

0.57258765x10“^ 

0.36964876x10“^ 

-0.32978152x10“!^ 

-0.40774761x10“® 

-0.26191716x10“® 

0.16568211x10“!^ 

0.21032277x10“® 

0.13454055x10“® 

-0.58782654x10“® 

-0.76069363x10“® 

-0.48484250x10“® 

0.13880285x10“!! 

0.18213056x10“!! 

0.11569505x10“!! 

-0.19485870x10“!! 

-0.25815323x10“!! 

-0.16344262x10“!! 

0.12253427x10“!® 

0.16330410x10“!® 

0.10302405x10“!® 
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TABLE XVII. Coupling constants Xij = Xji for the rank-3 ^P 2 potential. 

All =-0.41423805x101 
Ai 2 = 0.13294740x102 
Ai3 =-0.16229175x102 
A 22 =-0.46727996x102 
A 23 = 0.58280329x102 
A 33 =-0.73111775x102 
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FIG. 12. Phase shifts for ^P 2 state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XVIII. Parameters an,m of the form factor of Eq. (4.4) for the ^Di potential. The 
first column lists the rank-1 parameters with An = 0.38702146x10^. 



02,m 

03,m 

0.89936516 

0.26418769x101 

0.23149927x101 

-0.26186344x10^ 

-0.76699755x101 

-0.69187833x101 

0.40420417x10^ 

0.10171559x102 

0.80438855x101 

-0.37254532x10^ 

-0.77481279x101 

-0.53201948x101 

0.24045068x10^ 

0.40778272x101 

0.23638415x101 

-0.11674331x101 

-0.15957158x101 

-0.74115241 

0.44371340 

0.48775235 

0.16543078 

-0.13423928 

-0.11887084 

-0.23546429x10-1 

0.32471944x10-1 

0.23244826x10-1 

0.82565419x10-3 

-0.62682883x10-2 

-0.36286167x10-2 

0.57391815x10-3 

0.96078361x10-3 

0.44617253x10-3 

-0.17882613x10-3 

-0.11613293x10-3 

-0.42259918x10-^ 

0.30875386x10-^ 

0.10972859x10-^ 

0.29721861x10-3 

-0.36884776x10-3 

-0.80105503x10-^1 

-0.14441536x10-° 

0.32121255x10-° 

0.44458731x10-'^ 

0.39214969x10-3 

-0.20613138x10-^ 

-0.18322050x10-3 

0.15747854x10-1° 

0.96582033x10-° 

0.54076554x10-1° 

-0.61093792x10-11 

-0.32154572x10-1° 

-0.10765494x10-11 

0.24518321x10-12 

0.72109404x10-12 

0.12907316x10-13 

-0.45959801x10-1^ 

-0.97789252x10-1^ 

-0.70194793x10-1° 

0.35260306x10-1° 

0.60733413x10-1° 










TABLE XIX. Coupling constants = Xji for the rank-3 ^Di potential. 

All =-0.12403227x10® 

Ai 2 = 0.72374384x10^ 

Ai 3 =-0.33925562x10^ 

A 22 =-0.42322025x10^ 

A 23 = 0.19982358x10^ 

A 33 =-0.95274067x10® 
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FIG. 13. Phase shifts for ^Di state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XX. Parameters an,m of the form factor of Eq. (4.4) for the ^D 2 potential. The first 
column exhibits the rank-1 parameters with An = 0.59024361x10^ 



0-2,m 

(^3,m 

0.57066711 

0.16484822x101 

0.14099209x101 

-0.16583375 xlO^ 

-0.47749159x101 

-0.42072799x101 

0.25409617x10^ 

0.62850269x101 

0.48779676x101 

-0.23332252 xlO^ 

-0.47697116x101 

-0.32222764x101 

0.15018139x101 

0.25024490x101 

0.14316862x101 

-0.72790438 

-0.97809048 

-0.45159892 

0.27639429 

0.29955475 

0.10332838 

-0.83595442x10-1 

-0.73538742x10-1 

-0.161848780x10-1 

0.20227441x10-1 

0.14601065x10-1 

0.127725469x10-2 

-0.39078222x10-2 

-0.23406505x10-2 

0.124271741x10-2 

0.59974333x10-2 

0.30037269x10-2 

-0.610339620x10-^ 

-0.72616864x10-^ 

-0.30435345x10-^ 

0.111994241x10-1 

0.68758985x10-2 

0.23895692x10-2 

-0.133580683x10-2 

-0.50325849x10-® 

-0.14163274x10-® 

0.113449284x10-® 

0.28016277x10-1' 

0.60927632x10-® 

-0.701058725x10-® 

-0.11587213x10-® 

-0.17749233x10-° 

0.313672923x10-° 

0.34341671x10-1° 

0.29770498x10-11 

-0.991409110x10-11 

-0.68698690x10-12 

-0.11337993x10-12 

0.210314951x10-12 

0.82832450x10-1^ 

-0.48655006x10-12 

-0.269557305x10-11 

-0.45347221x10-1® 

0.63794448x10-11’ 

0.158651456x10-1® 
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TABLE XXL Coupling constants = Xji for the rank-3 ^D 2 potential. 

All = 0.28519285x10^ 

Ai2 =-0.17108157x10^ 

Ai 3 = 0.84546586x10^ 

A 22 = 0.10247694x10^ 

A 23 =-0.50401399x10^ 

A 33 = 0.24628297x10^ 
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FIG. 14. Phase shifts for ^D 2 state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XXII. Parameters an,m of the form factor of Eq. (4.4) for the potential. The first 
column gives the rank-1 parameters with An = 0.48683315x10^ 



0-2,m 

(^3,m 

0.70053005 

0.20372307x101 

0.17586107x101 

-0.20371819x10^ 

-0.59058643x101 

-0.52507453x101 

0.31301669x101 

0.77955339x101 

0.60941222x101 

-0.28784057x101 

-0.59247131x101 

-0.40278443x101 

0.18547114x101 

0.31124541x101 

0.17900376x101 

-0.89957466 

-0.12173452x101 

-0.563742149 

0.34172242 

0.37272087 

0.128054107 

-0.10337205 

-0.91321241x10-1 

-0.195064865x10-1 

0.25011910x10-1 

0.18051377x10-1 

0.129543369x10-2 

-0.48311038x10-2 

-0.28709573x10-2 

0.243205612x10-3 

0.74115882x10-3 

0.36379018x10-3 

-0.950679435x10-4 

-0.89691472x10-^ 

-0.36150992x10-^ 

0.169973636x10-4 

0.84868416x10-3 

0.27543762x10-3 

-0.203045225x10-3 

-0.62064686x10-3 

-0.15550511x10-3 

0.174498404x10-3 

0.34516816x10-'^ 

0.61282950x10-3 

-0.109761505x10-'^ 

-0.14259032x10-® 

-0.14661306x10-° 

0.502037641x10-° 

0.42202574x10-1° 

0.10222286x10-11 

-0.162758205x10-1° 

-0.84290188x10-12 

0.50688694x10-13 

0.355038528x10-12 

0.10144471x10-13 

-0.15433806x10-14 

-0.468513806x10-14 

-0.55417256x10-13 

0.14088408x10-13 

0.283746588x10-13 
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TABLE XXIII. Coupling constants Xij = Xji for the rank-3 potential. 

All = 0.17437597x10^ 

Ai 2 =-0. 10347752 xlO^ 

Ai 3 = 0.50147622x10^ 

A 22 = 0.61393167x10^ 

A 23 =-0.29733363 xlO^ 

A 33 = 0.14387562 xlO^ 
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E,„[MeV] 

FIG. 15. Phase shifts for state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XXIV. Parameters an,m of the form factor of Eq. (4.4) for the ^E 2 potential. The first 
column gives the rank-1 parameters with An =-0.15280632x10^ 



0-2,m 

0-3,m 

-0.10519807 

-0.10795104x101 

-0.11657860x101 

0.26681674 

0.22362235x101 

0.20871658x101 

-0.38458973 

-0.29306786x101 

-0.25641087x101 

0.33899788 

0.22258120x101 

0.17311373x101 

-0.20956819 

-0.11544199x101 

-0.759879540 

0.96873287x10-1 

0.42971181 

0.208280611 

-0.34663036x10-1 

-0.11609305 

-0.224855812x10-1 

0.97602881x10-2 

0.21859658x10-1 

-0.890945285x10-2 

-0.21786339x10-2 

-0.24236811x10-2 

0.551995844x10-2 

0.38647216x10-3 

-0.54141248x10-3 

-0.161738180x10-2 

-0.54465017x10-^ 

0.60393519x10-1 

0.316664790x10-3 

0.60817595x10-3 

-0.12963497x10-1 

-0.447135268x10-1 

-0.53523779x10-6 

0.16249029x10-3 

0.466785193x10-3 

0.36794342x10-1 

-0.13891113x10-6 

-0.362189145x10-6 

-0.19476964x10-8 

0.83947022x10-8 

0.207493943x10-1 

0.77622991x10-1° 

-0.35813659x10-° 

-0.862103751x10-° 

-0.22468227x10-11 

0.10514171x10-1° 

0.251419959x10-1° 

0.44459234x10-13 

-0.20104072x10-12 

-0.485787313x10-12 

-0.53628444x10-13 

0.22387162x10-11 

0.556075874x10-11 

0.29648478x10-11 

-0.10924776x10-16 

-0.284419168x10-16 
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TABLE XXV. Coupling constants Xij = Xji for the rank-3 ^F 2 potential. 

All =-0.22056858x10^ 

Ai 2 = 0.53156082x10^ 

Ai3 =-0.29476815x10^ 

A 22 =-0.13088707x10^ 

A 23 = 0.73245374x10^ 

A 33 =-0.41218725x10^ 
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E,„[MeV] 

FIG. 16. Phase shifts for ^F 2 state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 


TABLE XXVI. Parameters an,m of the form factor of Eq. (4.4) for the ^Es potential. The first 
column lists the rank-1 parameters with An =-0.15533123x10^ 



0-2,m 

0-3,m 

-0.10355351 

-0.10619702x101 

-0.11465411x101 

0.26271214 

0.22009693x101 

0.20540811x101 

-0.37869541 

-0.28848496x101 

-0.25239699x101 

0.33382672 

0.21914494x101 

0.17046664x101 

-0.20638383 

-0.11368554x101 

-0.748665750 

0.95407758x10-1 

0.42333402 

0.205500714 

-0.34141350x10-1 

-0.11445199 

-0.223823711x10-1 

0.96143050x10-2 

0.21585398x10-1 

-0.867190389x10-2 

-0.21462844x10-2 

-0.24058411x10-2 

0.540397503x10-2 

0.38078040x10-3 

-0.75437723x10-6 

-0.158522915x10-2 

-0.53669961x10-^ 

0.58646658x10-1 

0.310495776x10-3 

0.59937965x10-3 

-0.12652563x10-1 

-0.438486050x10-1 

-0.52756756x10-6 

0.15883035x10-3 

0.457758155x10-3 

0.36271728x10-1 

-0.13585744x10-6 

-0.355154343x10-6 

-0.19202569x10-8 

0.82111556x10-3 

0.203430516x10-1 

0.76537251x10-1° 

-0.35024827x10-° 

-0.845015123x10-° 

-0.22155859x10-11 

0.10278242x10-1° 

0.246353323x10-1° 

0.43844231x10-13 

-0.19638883x10-12 

-0.475781212x10-12 

-0.52889751x10-13 

0.21844895x10-11 

0.544288742x10-11 

0.29241713x10-11 

-0.10642051x10-16 

-0.278159322x10-16 
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TABLE XXVII. The coupling constants Xij = Xji for the rank-3 potential. 


All =-0.22415203x10^ 
Ai 2 = 0.54047498 xlO^ 
Ai3 =-0.29976913x10^ 
A 22 =-0.13315252x10^ 
A 23 = 0.74528575x10^ 
A 33 =-0.41950064x10^ 
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02,m 

03,m 

-0.12296582 

-0.12705383x103 

-0.13761116x103 

0.31102049 

0.26178431x103 

0.24457908x103 

-0.44801941 

-0.34259243x103 

-0.29980301x103 

0.39459099 

0.25961039x103 

0.20157805x103 

-0.24377244 

-0.13430375x103 

-0.87938384 

0.11259879 

0.49776587 

0.23706582 

-0.40253696x10"^ 

-0.13336950 

-0.22922059x10-3 

0.11322223x10-^ 

0.24641103x10-3 

-0.11752433x10-3 

-0.25240891x10-2 

-0.25611382x10-2 

0.68583670x10-2 

0.44711789x10-3 

-0.68617181x10-^ 

-0.19848209x10-2 

-0.62915174x10-^ 

0.81559745x10-^ 

0.38692325x10-° 

0.70141424x10-3 

-0.16636723x10-^ 

-0.54555554x10-4 

-0.61631030x10-3 

0.20532442x10-3 

0.56956135x10-° 

0.42302672x10-'^ 

-0.17453341x10-° 

-0.44240794x10-° 

-0.22361094x10-3 

0.10535983x10-7 

0.25394314x10-7 

0.89005384x10-^° 

-0.45039190x10-° 

-0.10580995x10-° 

-0.25735098x10-^3 

0.13286312x10-3° 

0.30978486x10-3° 

0.50877050x10-33 

-0.25608145x10-32 

-0.60169351x10-32 

-0.61321127x10-33 

0.28865188x10-34 

0.69356947x10-34 

0.33875686x10-37 

-0.14345838x10-3° 

-0.35807821x10-3° 
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TABLE XXIX. Coupling constants Xij = Xji for the rank-3 ^F 4 potential. 

All =-0.18792638 xlO^ 

Ai 2 = 0.45033170x10^ 

Ai3 =-0.24918509x10^ 

A 22 =-0.11024146x10^ 

A 23 = 0.61549429x10^ 

A 33 =-0.34552195x10^ 
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E,^[MeV] 

FIG. 18. Phase shifts for ^F 4 state without Coulomb effects. The crosses denote the 
RGM+OCM+LS result, the dashed line the rank-1 result, and the solid line the rank-3 result, 
respectively. 

In the higher energy region (above the break up threshold for ^He—i>d-|-p), absorption 
effects stemming from inelasticity should be considered. This requires the introduction of an 
imaginary component in the potential. We modified our potential to reproduce the reflection 
parameters, given by Yoshino et al. at TlLa6=19.48 MeV, using the anzatz 

1/C(LS) _ ^ A*[Ci{L • S) + C2{L • Sf + C^{L • 5)3]}i/Mocm^ 

where the parameters Gq, Gi, G 2 , G 3 are the same as in the 1 /OCM(ls) potential. And the 
complex factor A is chosen as follows 

A = l + ia{l-e-^^^e{Eth) (4.7) 

where A* is the complex conjugate of A, and a=21.0, 6=2.50 x 10“®, respectively. The break 
up threshold energy Eth is given by Eth = E — E^ + E^, E^ and E^ are the ^He and the 
deuteron binding energies respectively, while 6{x) is the step function. Using this complex 
potential, we calculated the phase shifts and found that they change, at most, by a few 
degrees in the region of 19.48 MeV. Therefore, we do not present these results in this paper. 


V. DISCUSSION AND CONCLUSIONS 

We constructed the p-^He effective potentials up to L=3 based on the RGM combined 
with the OCM technique which removes the PFS. These models do not include either an LS 
or a tensor force, which were not included in the nucleon-nucleon potential used to construct 
the nucleon-trinucleon potential. In the present work, we introduce the LS interaction 
phenomenologically at an intercluster level. This successfully describes the degeneracy in 
the spin-triplet (^Po,^Pi, ^P 2 ), (^Di, ^ 02 ,^ 03 ), and (^F 2 , ^F 3 ,^F 4 ) channels. The tensor 
force could be similarly included. However, these forces could be safely omitted as their 
influence on the data was shown, in the recent work of Yoshino et al. [|^, to be very small. 

In the higher energy region (above the break up threshold), absorption effects can be 
included via Eq. ( |4.(j| ). Since, however, our main concern is the construction of potentials 
at low energies, little attention was given to the construction of complex potential (see forth 
Eq. (P)). 
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The theoretical (RGM) phase shifts are well reproduced by our potentials. As far as 
the experimental phase shifts are concerned, it should be noted that the results obtained 
by various analyses are not in agreement to each other. However, our potentials fit the 
most recent phase shifts quite well. Furthermore, we compared our results with those 
of Ref. |IO, which are based on a Yamaguchi type separable potential htted to lower energy 


scattering data. The latter potentials reproduce well the differential cross sections at very 
low energies. 

Special attention was paid to effects of non-central forces in each of the negative parity 
states, 2“, 1“, 0“. We derived the energy spectrum by using the phase-shift data by Bel- 
tramin et al. (BFP) Ref. employing the three-dimensional Spline function interpolation 
method. Our result is compared with the BFP spectrum and with the GCM calculation Ref. 
^ as well as with the experimental data. The results are plotted in Fig. 0. The hrst and 


the second lowest levels of the BFP crossed each other. The GGM calculation reproduces 
the ordering of four negative parity states but the spectrum is totally shifted to a higher 
energy region, while our results are in good agreement with the experimental data in which 
the lowest level 2“ is due to the ^P 2 state, the second lowest 1“ level to ^Pi, the 0“ level to 
^Po, and the excited 1“ level to the ^Pi state. 



FIG. 19. ^Li resonance energy level above the p-^He threshold. The experimental data are 
taken from Ref. [22], Theor. BFP are the results of Ref. [10], Theor. GCM by Ref. [21], and 
Theor. Ours are the results of the present work. 

From the overall results, we may conclude that our potentials (RGM-I-OGM plus phe¬ 
nomenological LS force) reproduce the scattering and resonance state data of the p-^He 
system well. Concerning the OGM results, we note that they are essentially on-shell equiv¬ 
alent to the RGM one. However, we found that there are some differences, albeit small, in 
the phase shift for the P and D waves. 

Finally, our EST separable potentials can be used in few-cluster systems where p-^He is 
involved. One such reaction is ^He(d,p)^He, which can be treated as a three-body system, p- 
n-^He, within the Faddeev integral equations. Moreover, they can be used for interpretation 
purposes by constructing local equivalent interactions which can provide us a further insight 
into the interaction between light clusters, their characteristics concerning shape, range 
etc., as well as their bound and resonance structure. The understanding of the interaction 
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between light clusters will pave the way for a better treatment of few-cluster systems and 
their role in nuclear reactions and primordial nucleosynthesis. 
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APPENDIX A: RESONATING GROUP METHOD AND THE 
ORTHOGONALITY CONDITION MODEL 


In what follows we outline for convenience the ROM for the p-^He system. The total 
wave function for the p-^He system is dehned by 


T = A[0cl^(o-, r) x{Rn - J^cl)], 


(Al) 


where A is the antisymmetrization operator, ^(cr, r) the isospin function, and x{Rn — Rcl) 
the relative wave function between the clusters, while 0 cl is the internal wave function of 
the clusters which is dehned by superimposing two Gaussian functions 


1 3 

0CL = exp[--ai - RclY 

^ i=i 
1 3 

+ cexp[--a2Xl(u - ^cl) 

^ i=i 


(A2) 


Here Rql is the center of mass of the ^He cluster and r* are the coordinates of nucleons 
within the ^He cluster. The parameters ai, a 2 , and c are taken from Ref. [0, namely, 
= 0.25 fm“^, 02 = 0.71 fm“^, and c = 3.17. The Hamiltonian is given by 


H = Ho + Y1 

i>j 


(A3) 


where Hq is the total kinetic energy and Uj are the interactions between the (i, j) nucleons. 
The relative function x{Rn — Rcl) can be evaluated variationally by considering 


(JT|i7-U'|T) =0, 


(A4) 


where E' is the total energy of the whole system. In this model, the total energy is separated 
into the internal energy of the cluster Eint and the relative energy between the colliding 
clusters E, i.e. E' = E + Eint- Using Eq. ( |A2|) in Eq. ([A^) one can obtain the following 
integro-differential equation. 


|-V^ - Voir) + e) x(r) = ^ K(r, r')x(r') dr', 


(A5) 
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where Vd is the Direct local potential and K{v, r') is the corresponding nonlocal one which 
consists of three terms 

K{r, r') = Rxir, r') + Ky{r, r') + E'Mir, r'), (A6) 

where Kt is obtained from the kinetic energy term, Ky from the potential term, and M is 
often referred to as the norm-integral kernel or the norm kernel. In operator form Eq. ( |A5| ) 
is written as 


(ho + VD + KT + Ky + E'M)x = Ex (A7) 

where E and ho are the relative energy and the kinetic operator between the ^He cluster 
and the proton, respectively. The RGM (effective) potential can be identihed to 

^RGM(^) ^ =w + EM (A8) 

with 

w = VD + KT + Ky + EintM. (A 9 ) 

The normalization of the total wave function of the system is given by 

1 = (^1^) = MaMAMbx]) = {xl(i - AA)|x) = (V'lV'), (Aio) 

where 0a is the ^He cluster intrinsic wave function while 06 is a single nucleon wave function. 
It is, however, known that the norm of the relative function y is not one and therefore one 
dehnes the function 0, 


0 = Vl-Mx, (All) 

which has the proper probability density interpretation as the physical wave function. 

An additional problem with the RGM formalism is the existence of PFS. These states 
can be removed using the OGM technique in which Eq. (ATI) is rewritten as follows. 


(/!„ +Vb + Kt + Ky + E,„,M)x = £(1 - M)x- 


(A12) 


Furthermore, using Eq. (AHp one obtains 
1 


:(ho -|- Vd -|- Kt -|- Ky -|- EintM)- 


Therefore, the Pauli corrected intercluster potential is 


0 = E0. 


V 


OCM 


1 1 

{fiQ Vd + Kt + Ky -\- EintM ) — ,. . , — ht 


m-m 


01-70 


0 • 


(A13) 


(A14) 


A Pauli forbidden state is removed when the corresponding eigen-value of 01 — M vanishes. 
It should be noted that in this representation the energy dependence of the potential is 
eliminated. 
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The potential in Eq. (|A14|) is given in operator form. It can, however, be easily 

brought into a more suitable form for numerical calculation. For this, when the cluster wave 
function is given by one Gaussian term, the eigen-function of the norm kernel M is found 
using the harmonic oscillator function In the present case, however, where the ^He 

wave function is given by two Gaussian terms, the norm eigen-function is a superposition of 
harmonic oscillator functions. Now, the norm kernel M is given by 


where 


A/'= \Ut)G.j{U,\ 

ij=l 




(A16) 


(A16) 


The matrix elements of the norm kernel analytically while the matrix can be diagonalized 
numerically using the Jacobo method. For this, Eq. ( |A15|) is written as 


■^ = \Ui)aik'ykalj{Uj\ 

ijk=l 
oo 

= E \Ukhk{Uk\ 

k=l 


(A17) 


where 


\Ut) =T.\Ui)“it 


(A18) 


2=1 


Furthermore, the term 7 ^, is the eigen-value of the norm kernel and Uk is the corresponding 
eigen-function. Eq. (|A14|) then becomes 


where 


ij=l V-*- 


% 


■■{hoij + Wij) 


\/l - 


hoij \ {Uj 


hoij — (t/j I/iq I Gj ), 

Wij = {u.muj). 


(A19) 


(A20) 

(A21) 


APPENDIX B: MOMENTUM REPRESENTATION 


In what follows we also present the above potential in momentum space suitable for 
the Alt-Grassberger-Sandhas (AGS) equations . The momentum representation is intro¬ 
duced using the Fourier transforms. 


= 47 r / fP ■ je{kr) ■ dr , 

Jo 

= 47 r / ■ je{kr) ■ 3 i{k'r') dr dr'. 

Jo 


(Bl) 

(B2) 
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The potential consist of the following set of functions pO 
/j°^(a;r) = exp(—ar^) 


(B3) 

f^\a,a' : r, r') = (—exp(—ar^ — ar'‘^)je{ibrr') (B4) 

b : r, r') = (—i)Vr'(r^ + r'^) exp(—ar^ — ar'‘^)ji{ibrr') (B5) 

ff‘\a^b : r,r') = (—i)^rr'exp(—ar^ — ar'‘^){{ibrr')jLj^i{ibrr') — L ■ j(,{ibrr')} (B6) 

fi^\a,a' : r, r') = rr'exp(—ar^ — aV'^) (B7) 

fi^\a : r, r') = rr'(r^ + r'^) exp(—ar^ — ar'^) (B8) 

a', b, c, d : r, r') = {—iYrr^cr"^ + cV'^) exp(—ar^ — a'r'‘^)ji{ibrr') (B9) 


fP{a,d,b : r,r') = —brr'exp{—ar‘^ — a'r'‘^){i^'^^rr'jL+i{—ibrr') + i^^jei—ibrd)} (BIO) 


Using the dehnitions (|B1| ) and 

,71 


we obtain 


2a , + ,kk'. 

HfP) = 


-1 


^{fP} = TT^—UK^'^L+ii.a, a,b: k, k') 


^ a'k'^ + afc'^ 

e^P(-173- —)Jl+i/2{- 


bkk' 


Aaa' — 6^ 


4aa' — 6^' 


+ [ 


(4a^ — 6^)2 
4a(2L + 3) {Aa'^ + b'^){e + k'^) 


^{/f} = -|L + 


4a^ — (4a^ — b'^y 

(2L + 3)62 2abyk^ + k'^) 


Aa^ - 62 
(4a2 + b^bkk' 


(4a2 - 62)2 


CJjCJj 

^2 


2 1 


4a2 — 62)2 
X 2 ,+i(a, a,b : k, k') 
a'k"^ + ak'"^ , 


]X£(a, a,b ■. k, k') 
]X£(a, a,b k, k') 


Aaa' 


(5)i _ ro + ^'^)i 




4a 


, A :2 + A;' 2 ^ 
exp(-^-) 


4a 


^rU6)| ^ . r 4 a(^ + 3 / 2 ) _ 4a^2fc2 ^ ^2^/2 , 4 a(X + 3 / 2 ) _ 4a2/c^2 ^ ^2^2 . , . 

^ ^ ^ 4aa' - 62 (4aa' - 62)2 i + i 4^^/ _ 52 (4^^/ _ 52)2 iJ H , , • , ) 

Abkk'{a'c + ad) 


(4aa' — 62)2 


-lL+i{a,a',b : k, k') 


(Bi; 

(BU 


(Bi: 


(Bl^ 

(BU 

(Bl( 


(BU 


,26(X + 3/2) L 2b(a'k‘^ + ak'y, ^ . , , , , ,x kk'iAaa' + by , , , , , 

^ ^ 4aa'-62 + 6 “ ^ ^ ^ ^ ^ ^ ) • (^1^ 


(4aa' — 62)^ 


(4aa' — 62)" 


In the above 


Xya,a\b-.k,k') = J^{fYY 


(B19) 


while is the Hyperbolic Spherical Bessel Function related to the Spherical Bessel Function 
it by 


Ji+i/ 2 {x) = dxje{ix). 


(B20) 
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